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ABSTRACT 

We present a catalog covering 1.62 deg 2 of the COSMOS/UltraVISTA field with PSF-matched 
photometry in 30 photometric bands. The catalog covers the wavelength range 0.15/im — > 24/im 
including the available GALEX, Subaru, CFHT, VISTA and Spitzer data. Catalog sources have 
been selected from the DR1 UltraVISTA K s band imaging that reaches a depth of K a< t t — 23.4 AB 
(90% completeness). The PSF-matched catalog is generated using position-dependent PSFs ensuring 
accurate colors across the entire field. Also included is a catalog of photometric redshifts (z p hot) 
for all galaxies computed with the EAZY code. Comparison with spectroscopy from the zCOSMOS 
10k bright sample shows that up to z ~ 1.5 the z p hot are accurate to Az/(1 + z) = 0.013, with a 
catastrophic outlier fraction of only 1.6%. The z p hot also show good agreement with the z p hot from 
the NEWFIRM Medium Band Survey (NMBS) out to z ~ 3. A catalog of stellar masses and stellar 
population parameters for galaxies determined using the FAST spectral energy distribution fitting 
code is provided for all galaxies. Also included are rest-frame U-V and V-J colors, L 2 soo and Lir. 
The UVJ color-color diagram confirms that the galaxy bi- modality is well-established out to z ~ 2. 
Star- forming galaxies also obey a star forming "main sequence" out to z ~ 2.5, and this sequence 
evolves in a manner consistent with previous measurements. The COSMOS/UltraVISTA K s -selected 
catalog covers a unique parameter space in both depth, area, and multi-wavelength coverage and 
promises to be a useful tool for studying the growth of the galaxy population out to z ~ 3 - 4. 
Subject headings: infrared: galaxies - galaxies: evolution - galaxies: high-redshift - galaxies: funda- 
mental parameters 
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1. INTRODUCTION 

The last decade has seen significant progress in 
our understanding of the evolution of massive galaxies 
(Log(M star /M©) > 11.0) at < z < 5. These galax- 
ies exist in significant numbers already at 2 < z < 5, 
and their number density evolves fairly quickly over this 
redsh i ft ran g e (e.g.jMcLure et a l. 2006; Marchesin i et al.l 
[20091 l20lot iCirasuolo et al.l 120101 iCaouti et al.l \20lW . 
This early and rapid formation of massive galaxies 
is not well-reproduced in current theoretical models, 
where high-mass galaxies typically form much later (e.g., 
Fontanot et~aT1 120091: fGuo et al.l 120111 : iHenriques et all 
20H iBower et al.ll20rl . 

Between < z < 2 the growth of massive 

1 Based on data products from observations made with ESO 
Telescopes at the La Silla Paranal Observatory under ESO 
programme ID 179.A-2005 and on data products produced by 
TERAPIX and the Cambridge Astronomy Survey Unit on behalf 
of the UltraVISTA consortium. 

2 Catalog and other data products are available a t 
http: / /www. strw.leidenuniv.nl/galaxyevolution/ULTRAVISTA/ 
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galaxies is more gradual (e.g.. IBundv et al.l [ 2006: 
Arnouts et al.ll2007t iMarchesini et al.ll2009t lllbert et al.1 



2010HBrammer et all 1201 lUBielbv et al.ll2012t ). perhaps 
because this redshift range is also marked by a signifi- 
cant decrease in the star formation rates of these galax- 
ies. Current data show that the bi-modality between 
quiescent galaxies and star forming gal axies is largely 
estab l ished between 1 < z < 2 (e.g., [W illiam set al.1 
120091: IWhitaker et~aLll2011t iBrammer et al.ll2011l) . bT 
low z < 1, it appears the transformation of massive 
star forming galaxies into massive quiescent galaxies is 
largely completed and any subsequent mass growth in 
the massive population is primarily driven by mergers 
fe.g^ IBundv et alJl2006l : I A rnouts et al.ll2007tlijpert et al.1 
120101 IBrammer et alJl201lT ). 

This impressive recent progress in measuring the 
evolution of massive galaxies has primarily been 
made possible by the extensive investments made in 
ground-based near infrared (NIR) imaging capabili- 
ties over the last decade. The early deep pencil- 
beam NIR surveys that allowed for the fir st identifica- 
tions of massive, high-redshift galaxi e s (e.g..lFranx et all 
20031: iForster Schreiber et all 12004 iDaddi et all 120051: 



van Dokkum et_al.l 120061 ) such as FIRES (0.002 deg 2 . 



Labbe et al.l [2003h . G OODS (0.04 deg 2 iWuvts et all 
2008D and MUSYC (lOuadri et all [20071 0.12 deg 2 ), 
have now been superceded by much wider surveys 
that are equivalently deep, or de eper. These 



sur- 



vey^, such as the NMBS ( 0.44 deg 2 , Ivan Dokkum et all 



20091 IWhitaker 



et al.ll2011h. the UKIDSS-U DS (0.7 deg 2 
Lawrence et al.ll2007t: IWilliams et al.ll2009D. and WIRDS 



(2.1 deg 2 .lMcCrac ken et al.ll2010t iBielbv et al.ll2012l ) are 



2 



now the best-studied cosmic windows onto the massive 
galaxy population and have provided the source data for 
most of the results previously mentioned. 

The latest in this series of wider and deeper NIR 
sky s urveys is the UltraVISTA survey (jMcCracken et ail 
120121 1 . UltraVISTA has imaging in four broadband NIR 
filters (YJHK S ) as well as one narrow band filter centered 
on Ha at z = 0.8 (NB118). UltraVISTA is the deepest 
of the VISTA public surveys, and when fully complete 
will cover an area of 1.8 deg 2 down to K s ~ 24.0, with a 
deeper component (referred to as "ultra-deep") coverin g 
0.75 deg 2 down to K s ~ 25.6 (see lMcCracken et al.ll2012h . 
The current first data release of UltraVISTA is based on 
approximately on e season of observing ti me and is now 
publicly available (|McCracken et aL| [2012). In those data 
the imaging reaches a 5cr depth of K s < 23.9 AB in a 2" 
aperture. The full-depth UltraVISTA dataset will be ac- 
quired over a period of 5 - 7 years. 

In addition to an unique combination of area and 
depth, one of the main strengths of Ultra VISTA is that 
the survey field is the COSMOS field ([Scoville et al.1 
2007). COSMOS has arguably the most impressive array 
of multiwavelength coverage of any degree-scale part of 
the sky. It contains X -ray data from XMM — Newto n 
([Hasinger et al.ll2007H . and Chandra (lElvis et al.ll2009D 
UV i maging from the GALEX satellite (|Martin et al.l 
[20051) . extensive optical broad-band and optical medium- 
band imaging from the CFHT and Sub aru telescopes 



TABLE 1 

Summary of Photometric Data 



(|Taniguchi et al.l 120071 iCapak eta! 



120071). mid-infrared 



[20071: 



Fraver et al 



data from Spitzer ([Sanders et al 

2009) , submillimeter data from Herschel (| Oliver et al 
20121) . millimeter data from AzTEC and MAMBO 
(|Scott et all 120081 lAretxaga etal] 1201 ID .as well as ra- 
dio o bservations from the VLA ( Schinnerer et al.l 120071 

2010) . The field also has spectrosc opy for ~ 25,000 
galax ies from the zCOSMOS-bright (ILillv et alJ I I2QQ7L 
120091) and zCOSMOS-deep surveys (|Lillv et al.l I2007h . 
In addition to this, i t also has Hubble ACS data 
(|Koekemoer et al.l [2007) covering the full field, as well 
as deep WFC3 NIR imaging and grism spectroscopy 
in part of t he field from the CANDELS and 3D- 
HST surveys (iKoekemoer et al.ll201lUGrogin et al.ll2011t 
iBrammer et al.ll2012D .~ This extensive multi-wavelength 
coverage makes COSMOS an attractive field for perform- 
ing studies of distant galaxies. 

In this paper we present and make available a 30-band 
photometric catalog of the COSMOS field covering the 
wavelength range 0.15/im — »■ 24/im. The selection of 
sources has been made using the high image-quality Ul- 
traVISTA K s -band which allows for efficient selection of 
mass-complete samples of galaxies up to z — 4. In addi- 
tion to the photometric catalog, we also make available 
a catalog of photometric redshifts, stellar population pa- 
rameters, rest-frame colors, and UV/IR luminosities for 
all galaxies in the survey. The stellar mass function of 
galaxies to z = 4 determined using this catalog is pre- 
sented in a companion paper ([Muzzin et al.|[2013h . 

The layout of this paper is as follows. In § 2 we discuss 
the various observational datasets included in the photo- 
metric catalog. In § 3 we outline the steps of the catalog 
creation such as the PSF matching, source detection, and 
data quality controls. In § 4 we present the catalog of 
photometric redshifts and compare these to other pho- 
tometric and spectroscopic redshift measurements. In § 
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PSF stars in each of the 9 subfields. The depths within the 2.1 aperture are on 
the PSF-matchcd images and hence are effective depths for the purpose of color 
measurements in the catalog. The quoted depth at 24/im is a total flux. 

5 we present a catalog of stellar population parameters, 
and examine the basic properties of galaxies in the cat- 
alog such as the galaxy bi-modality and star formation 
"main sequence" . We conclude with a summary in § 6. 
Throughout this paper we assume an J7a = 0.7, ft m — 
0.3, and Ho = 70 km s _1 Mpc -1 cosmology. All magni- 
tudes are in the AB system. 

2. DESCRIPTION OF THE DATASET 
2.1. Photometric Datasets 

There is a wealth of imaging data at various wave- 
lengths available for the COSMOS field. For the K s - 
selected catalog we have chosen to include 30 photomet- 
ric bands that cover the wavelength range 0.15/im — > 
24/im. These bands and the papers that describe these 
datasets are summarized in Table 1. 

The catalog is ba sed on the YJHK ? NIR imaging data 
from UltraVISTA ([McCracken et alj|2012j) . and the in- 
clusion or exclusion of the available datasets was chosen 
to provide a match in areal coverage and depth of those 
data. The optical data consists of broadband data taken 
with Subaru/SuprimeCam (g + r + i + z + B ,Vj), as well as 
u* data from the CFH T/MegaCam (|Taniguchi et al.l 
120071 ICapak et all 120071) . We have also included the 
12 optical mediu m bands (IA427 - IA827) from Sub- 
aru/SuprimeCam ([Capak et al.| [2007) for a total of 23 op- 
tical/NIR bands. Obser vations from the GA LEX FUV 
and NUV channels (see iMartin et al1l2005[ ). as well as 
the 3.6/«ti, 4.5/im, 5.8/Ltm, 8.0/Ltm and 24/im channels 
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Fig. 1 — Coverage diagram for the UltraVISTA/COSMOS field. 
The filled gold region shows the area covered by the K s -selected 
catalog. The catalog contains objects only in regions where both 
NIR and optical coverage is available. It covers a total of 1.62 deg 2 , 
when regions contaminated by bright stars are excluded. 

from Spitzer's IRAC+MIPS cameras fsee lSanders et ail 
2007) have also been included using a source-fitting tech- 
nique designed for determining robust colors in highly- 
blended imaging data (see § 3.5). 

Several datasets have not been incorporated into the 
catalog. These include the deep u*g'r'i'z' imaging taken 
as part of the CFHTLS-Deep survey. These data are 
deeper than the Subaru broadband data; however, they 
cover ~ 0.9 deg 2 near the center of the COSMOS field, 
which is only slightly more than half of the UltraVISTA 
area. In principle, these data should improve the quality 
of the photometric redshifts and stellar mass measure- 
ments for galaxies in that region of the survey; however, 
including them would cause the uncertainties in those 
quantities to be location-dependent. This makes it much 
more difficult to estimate uncertainties in quantities de- 
rived from the entire survey such as the stellar mass 
function or the color distribution of galaxies. For the 
same reason we have chosen to not inclu de the deep NIR 
medi um band imaging from the NMBS (jWhitaker et al.l 
120111 ). which covers only 0.22 deg 2 of the field. 

We have also not i ncluded the J-band imaging of COS- 
MOS from KPNO (jCapak et al.l [20071), nor the H and 
K„ im aging from CFHT / WIRC AM (jMcCracken et al.l 
l2010f) . Those data have full coverage of the field; how- 
ever, they use similar filters as those in UltraVISTA and 
are shallower. 

2.2. Field Geometry 

The geometry of the COSMOS field is roughly a square 
patch on the sky; however, the coverage from the vari- 
ous datasets does not overlap perfectly. In Figure 1 we 
plot a schematic view of the layout of the datasets used 
in the photometric catalog. The selection band for the 
catalog is the UltraVISTA K s -band which is shown as 
the black outline in Figure 1. The UltraVISTA full field 
covers ~ 1.8 deg 2 ; however, as Figure 1 shows it is placed 



slightly to the west of the centroid of the Subaru broad 
and medium band imaging (orange outline). This offset 
in the UltraVISTA field was necessary to ensure that the 
four "strips" that are observed as part of the ultra-deep 
component of UltraVISTA did not coincide with the po- 
sitions of bright stars. The effective area of overlap be- 
tween the Subaru optical data and the UltraVISTA NIR 
data is 1.62 deg 2 once bright stars have been masked, 
and this region is shown as the gold shaded region in 
Figure 1. All sources in the catalog are contained within 
this overlap region. 

As Figure 1 shows, there is complete coverage for 
the 1.62 deg 2 catalog region from GALEX, IRAC, and 
MIPS. The majority of the ACS coverage is also within 
the UltraVISTA area. 

3. PHOTOMETRY AND SOURCE DETECTION 

Measuring the colors of galaxies accurately is 
paramount for determining properties such as photomet- 
ric redshifts and stellar masses. The available imaging 
data in the COSMOS/UltraVISTA field comes with a 
wide range of Point Spread Function (PSF) shapes and 
sizes and these variations need to be accounted for in the 
color measurements. The range of image qu alities for the 
vario us filters is listed in Table 1 (see also ICapak et al.l 
2007). As Table 1 shows, the image quality ranges from 
as good as 0.5" full width half maximum (FWHM) in 
some of the optical bands, to as poor as 4" - 5" FWHM in 
the GALEX and MIPS 24/im imaging. The optical and 
NIR imaging have PSFs that are comparable in shape 
(i.e., roughly Gaussian) and FWHM (0.5" - 1.2") and 
therefore PSF matching of those bands is performed in 
a similar way using standard techniques (see § 3.1). The 
space-based imaging from GALEX, IRAC and MIPS 
have more complicated PSF shapes that have consider- 
able wings, as well as a much larger FWHM. Photometry 
for those bands is performed separately using a source- 
fitting code designed to measure accurate photometry for 
highly-blended sources and is described in § 3.5. 

3.1. Optical and NIR PSF Matching 

PSF matching between the optical and NIR bands is 
performed by degrading the image quality of all bands to 
the image quality of the worst-seeing band. This process 
is performed separately in different regions of the sur- 
vey. Region- dependent PSF matching is needed because 
full coverage of the COSMOS field required 9 Su prime- 
Cam pointings per filter (see ICapak et alj 1200 7). The 
result of these multiple pointings is that in addition to 
the image-quality variations between filters, there are 
also image-quality variations from pointing-to-pointing 
within a given filter. The range of PSF FWHMs for each 
filter is listed in Table 1. 

To perform the PSF matching we divided the survey 
into 9 patches (labeled COSMOS- 1 - COSMOS-9) closely 
tied to the positions of the SuprimeCam pointings. We 
note that this approach does not account for intra-stack 
seeing variations in the NIR bands which are caused by 
the sparsely-filled nature of VIRCAM mosaic. The intra- 
stack PSF variation s of the order a few hund redths of an 
arcsecond (see e.g.. iMcCracken et al.l 120121) . and hence 
are much smaller than the PSF differences between dif- 
ferent bands. 

Within each region 10 bright, unsaturated reference 
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Fig. 2. — Top left: Curve of growth for bright PSF stars in the field. Subaru optical broad bands (dashed lines), NIR bands (solid 
lines) and Subaru optical medium bands (dot-dashed lines) are shown. Top right: Same curves of growth normalized relative to the worst 
seeing image (IA464, ~ 1.05"). The color aperture is shown as the dotted vertical line. Bottom panels: Same as top panels but after PSF 
homogenization has been performed. Before PSF homogenization the range of fluxes within the color aperture is a factor of 1.4. After PSF 
homogenization colors within the color aperture are accurate to ~ 1%. 



PSF stars were chosen and co-added into a reference PSF 
for that filter/region. Stars have a range of colors, and 
therefore it was difficult to find PSF stars that were suit- 
ably bright but unsaturated from the u* — > K s bands. 
This was dealt with by choosing one set of stars for the u* 
— > IA679 bands, and a separate set of stars for the IA738 
— > K s bands. An additional set of PSF stars was chosen 
for the Subaru i + band. The i + imaging has superior im- 
age quality (~ 0.5") across the field, and was taken with 
longer exposure times than the other band s so all stars 
down to i + < 21.8 are saturated (see also ICapak et al.l 
2007). Finding unsaturated PSF stars required going to 
much fainter stars than in the other bands. 

With a reference PSF constructed for each band/region 
we used the task lucy from the IRAF STSDAS package 
to compute the convolution kernels necessary to degrade 
the images to the image quality of the worst-seeing band. 
In 2/9 of the regions the worst image-quality band was 
the IA464 band, with typical seeing of ~ 1.05". In the 
other 7/9 regions the worst image-quality band was the 
g + band with seeing of 1.0" - 1.2". 

In Figure 2 we plot a illustration of the PSF match- 
ing process for the COSMOS-1 field. In the upper left 
panel we plot the growth curve of the PSF stars before 
PSF homogenization is performed. In the upper right 
panel we plot the same growth curves normalized rela- 
tive to the worst-seeing band (IA464) . The dotted verti- 
cal line shows the 2.1" aperture used for color measure- 



ments. Prior to PSF homogenization, the dispersion in 
the flux contained within the color aperture for the best 
and worst image-quality bands is a factor of ~ 1.4. 

In the bottom panels of Figure 2 we plot the reference 
stars after PSF homogenization has been performed. As 
the right panel shows, the dispersion in flux within the 
color aperture for the reference PSFs is reduced to < 1%. 

3.2. Source Extraction and Photometry 

Source detection and photometry is performe d using 
the SExtractor package (|Bertin fc Arnoutslll996f ) in dual 
image mode with the non-PSF-matched UltraVISTA K s - 
band image as the detection image. Objects are detected 
by flagging pixels > 1.7c above the background, after a 
2x2 pixel convolution kernel has been applied. We re- 
quire 10 connected pixels for a detection. At face-value 
this appears to be a conservative detection limit; how- 
ever, it is required to prevent the detection of a signifi- 
cant number of spurious sources. The native pixel scale 
of VIRCAM is 0.34" pixel" 1 but the data have been 
resampled to the COSMOS optical data pixel scale of 
0.15" pixel" 1 (see iMcCracken et al.lf2012h . This makes 
the noise in the UltraVISTA images highly correlated 
causing SExtractor to underestimate it and return a high 
fraction of spurious detections if a more relaxed object 
detection criterion is used. 

Photometry is performed in each of the 9 regions inde- 
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Fig. 3. — Left panel: The color - color space used to perform star/galaxy separation. One-third of objects with K s < 21.5 are plotted. 
Objects that are classified as stars using SExtractor (class.star > 0.95) are plotted in blue and objects classified as galaxies (class_star 
< 0.95) are plotted in red. There is a clear separation of stars and galaxies in the u* - J vs. J - K s color - color space. In the catalog, 
objects that lie below the lines are classified as stars, and those above are classified as galaxies. Right panel: Fraction of objects that are 
stars as a function of K s -band magnitude. 



pendently and afterwards these are merged into a master 
catalog for the entire field. For each filter we determine 
two fluxes, a Kron-aperture flux based on SExtractor's 
f lux_auto parameter, as well as the flux within a circu- 
lar color aperture. The color aperture was chosen to be 
f4 pixels in diameter, corresponding to 2.1" on the sky. 
This size makes the color aperture ~ 2x the FWHM of 
the worst-seeing image in each region. 

For each galaxy we supply a total K s -band magnitude 
(K Sj j 0t ) based on the f lux_auto parameter from SEx- 
tractor. This flux is t he flux meas ured within 2.5 times 
the Kron radius (Rk, [Kronl Il980f ) . The measured flux 
within 2.5 x should acco unt for > 96% of the total 
flux of the galaxy (|Kronll980h . We correct the f lux_auto 
to a total flux by measuring the growth curve of bright 
stars out to a radius of 8". Depending on the Rk , this 
correction is of the order 2 - 4%, as expected. 

In the catalog we remove all objects with magnitudes 
< 3c limit in the K s -band (K s = 24.35, 2.1" aperture), 
which results in a final photometric catalog of 262 615 
sources. 



for each filter, and we list those depths and the range of 
depths in Table 1. The quoted depths in Table 1 are the 
5a depth in the 2.1" color aperture, and hence are the 
effective depth of the photometry. Given the uncertainty 
in the absolute IRAC zcropoint, we have also included an 
additional 5% systematic error to the IRAC magnitudes 
in quadrature. 

Errors in K s , tot are calculat ed using the method de- 
veloped by lLabbe et al.l (120031) . which as also u sed by 
iQuadri et al.1 (|2007t ) and IWhitaker et al.l (|201lD . The 
background noise in an aperture of diameter N pixels 
will scale as N 2 in the limiting case of perfect correlation 
between pixels, and will scale as N in the limiting case 
of no correlation between pixels. We parameterize the 
background noise in an aperture as o~ ap = <r 1 aN /3 , where 
a i is the size of the aperture in arcseconds, and 1 < /3 < 
2. We perform the empty-aperture measurement for a 
range of aperture sizes and fit this relation for a and /3, 
finding f3 — 1.82 and a — -0.35. From this relation we 
then compute the error in the f lux_auto based on the 
rms noise in an aperture of 2.5R_r-. 



3.2.1. Galactic Extinction Correction 

The photometry in all bands is correct ed for Galactic 
dust a ttenuation using the dust maps from Schlcgcl et al. 
(1998). Dust extinctions are calculated in each of the 
9 photometry regions individually. The COSMOS field 
does not have substantial Galactic dust, and the correc- 
tions are of order 15% in the GALEX bands, 5% in the 
optical and < 1% in the NIR and MIR. 

3.2.2. Photometric Errors 

All of the optical and NIR images have been resampled 
during the mosaicking process. The resampling causes 
correlations in the noise between pixels and therefore 
standard Poission estimation of the background noise is 
not a reliable estimate of the photometric errors. We 
measure the background noise in each band, and in each 
region, by placing 10 000 empty apertures across the field 
and measuring the rms flux in those empty apertures. 
These rms estimates also give a measure of the 5cr depth 



3.3. Star Galaxy Separation 

Similar to the NMBS, star galaxy separation is per- 
formed use the J - K s vs. u* - J color space. In Figure 
3 we plot this color space for a randomly-selected sub- 
sample of 30% of the objects with K s < 21.5. Points 
are color-coded by SExtractor's class_star parameter, 
with star-like profiles in blue, and galaxy-like profiles in 
red. As Figure 3 shows, the two types of objects clearly 
segregate in color space. Based on this segregation, ob- 
jects are classified as galaxies if they meet the following 
color criteria, 

J -K s > 0.18 x (u* - J) - 0.75, [u - J < 3.0], (1) 

J — K s > 0.08 x (u* - J) - 0.45, [u — J > 3.0] (2) 

As shown in the right panel of Figure 3, stars dominate 
the catalog at K s j t < 17, but make up < 10% of objects 
at K Mot > 21.0. ' 
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3.4. Survey Completeness 

In order to construct mass- or luminosity-limited sam- 
ples, the completeness as a function of the selection band 
needs to be quantified. To estimate the point-source 
completeness we use the PSF stars as template sources, 
and insert these into the original non-PSF-matched K s - 
band image. We then attempt to recover these sources 
with SExtractor. When recovering the simulated stars 
we used the identical SExtractor parameters as were used 
for object detection in the catalog. 

This completeness test is performed in two ways. 
Firstly, we insert sources into a version of the image 
where all objects have been masked out using SExtrac- 
tor's segmentation map, and then attempt to recover 
these objects with SExtractor. The recovery rate of these 
objects gives an estimate of the completeness based solely 
on the noise properties of the image. We then perform 
the test again, this time using the real image with all 
objects still included. The recovery rate of this method 
gives a measure of the overall completeness, which is al- 
ways less than unity because of real effects that cause 
objects to be missed such as contamination from bright 
stars, or blending with other galaxies. 

In Figure 4 we plot the completeness curves as a func- 
tion of K Si tot- The 90% completeness limit in terms of 
the noise properties alone is K Sjtot = 24.0 mag; however, 
this corresponds to < 80% completeness in the real data. 
Based on the noise properties, the 100% completeness 
level of the UltraVISTA K s , tot data is 23.4 mag. This 
K Sj t ot limit corresponds to a 90% completeness level in 
the actual data. At magnitudes fainter than 23.4 the 
data becomes incomplete quite rapidly. This can be seen 
as the inflection point in the number counts i n the inset 
of Figure 4 (see also iMcCracken et all l2012f) , which is 
the result of the increasing incompleteness as well as the 
increase in photometric uncertainties due to poor S/N 
near the completeness limit. In order to construct mass- 
complete samples we only use sources in the catalog with 
K s ,tot < 23.4. We note that this completeness is 0.08 
mag brighter than wo uld be measured using th e DR1 Ul- 
traVISTA data from IMcCracken et all (|2012D . We have 
adjusted the zeropoint of the DR1 images by this amount 
in order to bring it into better agreement with other K s 
surveys such as 2MASS and the NMBS (see Appendix). 

Also shown in the inset of Figure 4 is t he sample com- 
plete ness measurement from the NMBS (Whi taker et al.l 
120111) . The 90% completeness of the NMBS is K tot < 
22.8, showing that UltraVISTA is ~ 0.6 mag deeper than 
NMBS. The surveys have similar exposure times on 4m- 
class telescopes, so the extra depth is primarily a result of 
the su perior image quality in the UltraVISTA DR1 [r 
0.75", IMcCracken et al.l 
(~ 1-1", 



Whitaker et al 



2012) compared to the NMBS 



20ll 



3.5. IRA C and MIPS Photometry 

The PSF FWHM of IRAC and MIPS are approxi- 
mately 2 - 5x larger than the FWHM of the worst- 
seeing optical and NIR PSF FWHMs. This means that 
degrading the image quality of the optical/NIR data 
to that image quality would cause a substantial reduc- 
tion in the signal-to- noise (S/N) of the color measure- 
ments as well as unnecessary blending of well-resolved 
galaxies. Therefore, instead of degrading PSFs we use 




19 20 21 22 23 24 25 
Ktot [AB] 

Fig. 4. — Completeness fraction as a function of total K s -band 
flux for point sources. The completeness is calculated using sim- 
ulations which insert point sources into the K s -band image and 
attempt to recover them. The grey curve is calculated by masking 
all real sources in the image and demonstrates the completeness 
based on the background noise in the image. The black curve is 
calculated by keeping in all real sources and therefore is lower be- 
cause of effects such as blending and confusion. In the inset we 
show the differential number counts in K s -band and a compar i- 
son to the number counts from the NMBS (Whitaker et al. 2011). 
Based on the simulations and the inflection in the number counts, 
the 90% completeness limit of the catalog is K Sj t t < 23.4. 

a source-fitting code developed to measure deblended 
photometry of heavily-confused images. This code is 
well-tested and has been used in many previous K s 



select e d catalogs (see lLabbe et al 



2007; March esini et a 
Whitaker et al 



120091: 



many p 

1 1200a n 



Wuvts et all 
iWilliams etal.1 120091 : 
Marchesini et al.l 1201 2f). includ- 
ing several with ultra-d eep IRAC imaging ()Labbe et al.l 
120101 lLabbe eTaI]l20ll) . 

In brief, we assume that there are no color gradients in 
galaxies between the K s band and the IRAC and MIPS 
bands. The K s band is then used as a high- resolution 
template image to deblend the IRAC and MIPS pho- 
tometry. Each source extracted from the K s image is 
convolved with a kernel derived from bright PSF stars in 
the K s and IRAC/MIPS images. The convolved galaxies 
are then fit as templates in the IRAC and MIPS bands 
with the total flux left as a free parameter. In this pro- 
cess, all objects in the image are fit simultaneously. Once 
the template fitting is converged, a "cleaned" image is 
produced for each object in the catalog by subtracting 
off all nearby sour ces (for an example of this process see 
iWuvts et al.ll2007l Fi gure 1). Aperture photometry is 
then performed on the cleaned image of each source. For 
the IRAC (MIPS) channels the photometry is performed 
in a 3" (5") diameter aperture for each object. This flux 
is then corrected to a flux within the 2.1" color aperture 



"IRAC 



(2.1' 



F, 



I RAC '.cleanei 



>(3'' 



K, 



conv, worst — PS F 



(2.1") 



K, 



conv ,1 RAC 



0") 



(3) 

where K C onv,worst-PSF (2-1") is the flux in the 
PSF-matched K s image in the color aperture and 
,irac(3") is the flux in the K s image convolved 
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to the IRAC PSF. Given that MIPS fluxes are primar- 
ily used to indicator SFRs, and are not used for colors, 
within the catalog the MIPS fluxes are listed as total 
fluxes rather than color fluxes. The MIPS fluxes have 
been converted to total fluxes using an aperture correc- 
tion of a factor of 3.7, as listed in the MIPS instrument 
handbook 10 . 

3.6. GALEX Photometry 

The image quality of the GALEX NUV and FUV data 
are also substantially poorer than that available from the 
ground-based data (~ 4" - 5"). We cleaned and photome- 
tered the GALEX data using the same source-fitting 
code as for the IRAC and MIPS data. For the GALEX 
data we use the CFHT u*-band data as a template rather 
than the K s -band because the u* band and the FUV and 
NUV bands have a close correspondence in wavelength. 
This required creating a new source list based on the 
u* images and then matching that source list to the K s 
source list. When performing this matching we use a 
conservative matching tolerance (< 0.5") to ensure only 
true counterparts are matched to the K s -selected cata- 
log. No GALEX photometry is provided for K s -selected 
sources that do not have a clear match to the u*-selected 
sources. We note that only the GALEX photometry is 
provided by matching sources, the u*-band photometry 
itself is still measured using the standard PSF matching 
method and SExtractor in dual-image mode. 

3.7. zCOSMOS Spectroscopic Redshifts 

We match the photometric catalog to the catalog of 
spectroscopic redshifts (z svec ) available fro m the zCOS- 
MOS lOk-bright sample (|Lillv et all 120091) . That cata- 
log contains ~ 10 000 redshifts for galaxies with i' < 
22.5. Due to the bright limit, the majority of redshifts 
in the catalog are at 0.1 < z < 1.5. The matching is 
performed for the subset of zCOSMOS sources with red- 
shift confidence classes in the range 3 < CC < 5. These 
are the highest-confidenc e spectroscopic r edshifts, and 
should be 99% accurate ([Lilly et al.ll2007l ). We match 
only the high-confidence z spec because these are used to 
determine offsets in the zeropoints in various filters (§ 
4.1), and to determine quantities such as stellar masses 
and rest-frame colors. Uncertain z spec can be extremely 
wrong resulting in unphysical results for this process. 

Matching is done using an 0.5" search radius. The 
optical i + -band and the UltraVISTA K s -band are well- 
registered so unique matches are found for all sources. 
Overall, there are 5105 z spec from zCOSMOS matched 
to sources in the K s -selected catalog. 

3.8. Corrections for Bright Stars and Bad Regions 

The COSMOS field is large enough that it contains 
regions with very bright stars. These bright stars create 
reflections and large diffraction spikes that make the pho- 
tometry for galaxies in those regions unreliable. In the 
catalog we provide a parameter, contamination which 
indicates whether an object's photometry has been con- 
taminated by a nearby bright star. The contamination is 
determined by first generating a source list of optically- 

10 http://irsa.ipac.caltech.edu/data/SPITZER/ 
docs /mips/ mipsinstrumenthandbook 



bright stars within the COSMOS field from the USNO- 
B catalog, as we ll as a list of NIR-brig ht stars from the 
2MASS catalog (|Skrutskie et all 12009 ). We determined 
an empirical relation between the brightness of the star 
and radius of the contaminated photometry All objects 
that are within the contamination radius of a bright star 
in either the optical or NIR are then given a contamina- 
tion flag = 1. 

In general we have been conservative with the contami- 
nation radius in order to provide the most reliable photo- 
metric catalogs possible. Still, the total area considered 
to be contaminated by bright stars is fairly modest. The 
full coverage region of the catalog covers an area of 1.68 
deg 2 , whereas the usable, uncontaminated region covers 
1.62 deg 2 . This means only ~ 4% of the total area is lost 
to bright stars. 

The Subaru optical bands also contain some small 
square-like regions throughout the survey that do not 
have data. Some of these regions are located at the cen- 
tral few pixels of bright stars, but others occur through- 
out the survey area, seemingly at random locations. It is 
unknown what the source of these regions are, but they 
occur most frequently in the i + band data. The pixel 
values in these regions have either been set to "nan" or 
10~ 31 , and failure to mask these regions causes nonsen- 
sical fluxes to be measured for nearby galaxies. The re- 
gions occur frequently enough that we developed a pro- 
cedure for identifying them. 

A pixel map of each optical band showing the location 
of the "nan" region is made. This map is then smoothed 
with a kernel that grows the size of the contaminated 
regions by a factor of ~ 2. When the final catalog is con- 
structed these mask regions are then checked for each 
object, in each filter. If the object lies within a contami- 
nated region for that filter, the flux in that filter is set to 
-99, and the parameter nan_contam is incremented by 1. 
This setting for the flux makes the photometric redshift 
fitting code (§ 4) and stellar population fitting code (§ 5) 
effectively consider the object not-observed in that par- 
ticular filter. This is useful because the nan regions are 
in different locations in various filters. The filter-by-filtcr 
masking allows us to effectively keep objects with some 
contamination in the catalog, but remove only the pho- 
tometry in the contaminated filters. Because nan_contam 
is incremented, it is a metric of the number of filters 
that have contamination for a given object. It is recom- 
mended to only use objects with contamination in < 5 
filters. 

3.9. Photometric Catalog Layout and Useage 

The layout of the full photometric catalog is summa- 
rized in Table 2. The catalog is presented as a set of 
fluxes in the 2.1" color aperture with an AB zeropoint 
of 25.0. These can be converted to AB magnitudes via 
m x = -2.51ogio(fx) + 25.0, where x denotes a given filter. 
Also given is the total K s -band magnitude from SExtrac- 
tor 's f lux_auto, which has been corrected to a total flux 
using the growth curve of the PSF stars. The flux in any 
filter can be converted to a total flux via, 

fx.tot — fx x _ „ — • (4) 

JKs 

The catalog contains the star indicator (= 1 for stars, 
= for galaxies), as well as the contamination and 



TABLE 2 
Summary of Photometric Catalog 



Column 


Parameter Name 


Description 


1 


id 


Object identifier number 


2, 3 


ra, dec 


Right Ascension and Declination in J2000 decimal degrees 


3, 4 


xpix, ypix 


Pixel position of object in the K 3 image 


5 


Ks.tot 


Total K s -band flux with additional aperture correction applied 


6 


cKs_tot 


Error m total K s -band magnitude determined from scaled empty apertures 


7-67 


X, eX 


Flux and error in filter X measured in a 2.1" aperture from PSF-matched images 


68 


K.flag 


SExtractor's FLAG output for the K s -band image 


69 


K_star 


SExtractor's CLASS-STAR output from the K s band image 


70 


K_Kron 


Kron radius in the K s -band 


71 


apcor 


Aperture correction that has been applied to FLUX_AUTO to determine Ks_tot 


72, 73, 74 


z_spec, z_spec_cc, z_spec_id 


Spectroscopic redshift, spectroscopic redshift quality flag, and ID number from the zCOSMOS catalog 


75 


star 


Star/galaxy indicator determined from color-color plot (star = 1, galaxy = 0) 


76 


contamination 


Indicates proximity to a bright star in the optical or NIR (contaminated = 1, uncontaminated = 0) 


77 


nan_contam 


Number of filters where object lies near gaps in photometry or near a saturated star 


78, 79 


orig_cat_id, orig_cat_field 


ID and tile of object in the tile catalog 


80 


USE 


Indicates galaxies with uncontaminated photometry and S/N > 5 



nan_contam parameters described in §3.8. Lastly, we in- 
clude a parameter USE. Selecting objects with USE = 1 
in the catalog selects the subset of objects with star = 
0, contamination = 0, nan_contam < 5, and K s < 23.9. 
The latter is the 5a depth of the survey in the color aper- 
ture. Objects with USE = 1 are considered to be galaxies 
with uncontaminated photometry with fluxes sufficiently 
bright that the photometry is still accurate. 

4. PHOTOMETRIC REDSHIFTS 

Photometric redshifts (z p hot) are calculated for a ll 
galaxies using the EAZY software (jl rammer "et~all2 008h 
EAZY determines the z p h ot for galaxies by fitting 
their SEDs to linear combinations of seven templates, 
six of which are derived from the PEGASE models 
(|Fioc fc Rocca-Volmerangel I1999D . as well as an add i- 
tional red template from the models of lMarastonl (2005). 
A detailed d escription of EAZY's fitting process can 
be found in iBrammer et all (|2008l ). EAZY also ac- 
counts for potential mismatches between the data and 
the templates via th e "template error function" (see 
IBrammer et all [2008). The template error function 
weights photometric points in the template fitting based 
on their implied rest-frame wavelength. In particular, 
when using the default template error function, measure- 
ments corresponding to the rest-frame NIR are down- 
weighted compared to the rest-frame optical given the 
curre nt uncertainties in models at th ese wavelengths 
(e.g., iMarastonl [20051 : iKriek et aLll2010D . 

Photometric redshifts were determined with EAZY pri- 
marily using the default set of parameters, although sev- 
eral optimizations were added after examination of the 
initial output. For the first run we used the 7 default 
templates, the vl.O template error function, the K. s ^ot 
magnitude prior, and allowed photometric redshift solu- 
tions in the range < z < 6. 

Comparison of the resulting z p hoto with the z spec 's from 
zCOSMOS was good, although there were some signif- 
icant outliers (~ 5%). Some experimentation showed 
that catastrophic outliers could largely be eliminated by 
running EAZY with a 5% systematic error included for 
all photometric bands. The 5% systematic helps reduce 
catastrophic outliers because the optical medium bands 
have narrow bandpasses and therefore problems in the 
photometry in consecutive bands can create very sharp 



features in the SED. 

After eye-examination of the best-fit models for a sub- 
sample of the galaxies it became apparent that there were 
two populations of galaxies that were not well-described 
with the default template set. The first of these were 
galaxies at z > 1 with post starburst-like SEDs. These 
galaxies tend to have strong Balmer breaks, but also a 
continuum that is very blue redward of the Balmer break, 
and it is not possible to reproduce such an SED with 
the default EAZY template set . We added a 1 Gyr old 
single-burst iBruzual fc Charlotl ()2003[ ) model to the tem- 
plate set and this significantly improved the fits for these 
galaxies. The ne ed for a similar temp late in EAZY was 
also reported bv lOnodera et all (|2012t ). 

The second set of galaxies with problematic fits was a 
population of UV-bright galaxies at 1.5 < z < 3.5 with 
UV continua that were blue, but still redder than could 
be produced with the default templates. Given that the 
catalog is K s -selected, these are likely to be the most 
massive part of the Lyman Break Galaxy (LBG) pop- 
ulation, and they may have slightly more dust extinc- 
tion than lower-mass LBGs. We added a slightly dust- 
reddened young population to the template set and this 
improved the fit to this population. 

4.1. Zeropoint Offsets 

Photometric redshifts are extremely sensitive to errors 
in photometric zeropoints. This is especially pronounced 
when medium bandwidth filters are used as they tend 
to be closely spaced in wavelength and errant zeropoints 
can create sharp features in the SEDs. A common pro- 
cedure to ensure the best-quality z p h oto is to refine the 
photometric zeropoints of a catalog using a subsample 
of ga l axies with spectroscopic redshifts (e.g-.lllbert et all 
I21MI2001 IBrammer et all[201ll iWhitaker et alJl20lTIT 7 

We perform this process using an ite rative zeropoint 
offset code developed for the NMBS (see Whitakc r et all 
120111 ). We use only the highest-quality spectroscopic 
redshifts from zCOSMOS (3 < CC < 5) as well as a 
set of 19 z spec of massive galaxies at z > 1 from other 
programs in the COSMOS field (lOnodera et all [2011 
Ivan de Sande et al.ll201lL l2012t iBezanson et alll2013[ ). 

The procedure is as follows. Firstly, EAZY is run on 
SEDs of the spectroscopic redshift sample fixing the pho- 
tometric redshift to the spectroscopic redshift. The me- 
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TABLE 3 
Zeropoint Offsets 
Determined from 
Spectroscopic Redshifts 



Filter 


Offset 




('Map-') 


fuv 


— 


nuv 


- 


u* 


-0.10 


B J 


-0.12 


g+ 


-0.11 


V, 
r+ 


0.03 


-0.06 


i+ 


-0.03 


z+ 


-0.06 


IA427 


-0.12 


IA464 


-0.12 


IA484 


-0.08 


IA505 


-0.07 


IA527 


-0.09 


IA574 


-0.17 


IA624 


-0.09 


IA679 


0.06 


IA709 


-0.10 


IA738 


-0.15 


IA767 


-0.13 


IA827 


-0.17 


Y 


-0.10 


J 


-0.14 


H 


-0.18 


K s 


-0.08" 


3.6/^m 




4.5/^m 




5.8/im 




8.0^m 




24/im 





Note. — Offsets are de- 
fined such that ZPeazy — 

^nominal + offset. 



"The K a zeropoint offset 
is determined by comparing 
photometry with other sur- 
veys, not from the spectro- 
scopic redshifts 

dian offset in each filter compared to the best-fit tem- 
plate is measured for the full sample. The photometry 
for all galaxies is adjusted by this amount and EAZY is 
re-run. The process is repeated until the median offsets 
converge to a value of less than 0.01 mag in every fil- 
ter. When calculating the offsets, the K s -band is used as 
the "anchor" filter and is not adjusted. Given that the 
offsetting is based on colors, having a filter that is not ad- 
justed is important so that the zeropoints do not drift in 
an absolute sense while the relative zeropoints are being 
improved. We note though that the K s zeropoint was ini- 
tially adjusted by 0.08 mag to bring it into better agree- 
ment with other surveys (see Appendix). We have also 
not computed offsets for the IRAC and GALEX bands. 
It is unclear how well the EAZY templates should re- 
produce these wavelength ranges and so we chose to not 
iterate them as it may introduce incorrect colors because 
of the adopted template set. 

In Table 3 we list the zeropoint offsets determined 
for all bands. These are typically small for the optical 
broad band filters (~ 0.05 mag), but can be larger for 
the optical medium bands (up to 0.17 mag). Similar size 
offsets for those band s were also found in the NMBS (see 
IWhitaker et~aT]|2011l Table 4). Offsets of the order 0.1 - 



0.2 mag are seen for the UltraVISTA bands. These are 
comparable to the offsets seen betwee n the UltraVISTA 
data a nd the CFHT WIRCAM data bv lMcCracken et all 
(l20ll . 

4.2. Photometric Redshift Catalog and Comparison with 
Spectroscopic Redshifts 

In Table 4 we show the layout of the photometric red- 
shift catalog from EAZY. The parameter z_peak corre- 
sponds to the peak probability of the P(z) function, and 
is considered to be the most likely z p hot- The 68% and 
95% confidence intervals are calculated by integrating the 
P(z) function. 

In the left panel of Figure 5 we show a compari- 
son between the z p hot and z spec for the zCOSMOS red- 
shifts , as well as the high-redshift spectrosc o pic sa mples 
from lOnodera et alj (l20ll) lBezanson et all ([2013D . and 
Ivan de Sande et al.1 ()2012f ). In general, the photometric 
redshifts compare well to the spectroscopic redshifts, al- 
though we note that this comparison is primarily from 
galaxies at z spec < 1.5. Galaxies that are > 3cx outliers 
from the one-to-one relation based on their redshift er- 
rors from EAZY are considered as catastrophic outliers. 
The fraction of these galaxies is low (1.56%), and the rms 
dispersion around the one-to-one relation for the remain- 
der of the sample is 0.013 in Sz/(l+z). 

The Zphot accuracy compares well to the z^ hot accuracy 
for other catalogs in the COSMOS field. Illbert et all 
(2009) measure a catastrophic outlier fraction of 0.7%, 
and 5z/(l+z) = 0.007. This is sl ightly bette r than , 
for our K s -selected catalog; however. Illbert et al.1 (2009) 
use a more restricted set of z spec for this test (~ 
4000 z svec compared to ~ 5000 z spec in our catalog). 
IWhitaker et al.l (|2011h find a catastrophic outlier fraction 
of 5% with a Sz/(l+z) = 0.008. This is a higher outlier 
fraction but a better rms dispersion. This difference is 
mostly likely from our inclusion of a 5% systematic error 
in the error bars when fitting with EAZY which reduces 
catastrophic outliers at the expense of some redshift pre- 
cision. 

4.3. Comparison with NMBS Photometric Redshifts 

The zCOSMOS z spec provide a useful diagnostic of the 
accuracy of the z p hot\ however, there are few z spec at z > 
1.5, leaving the accuracy in this redshift range less cer- 
tain. The NMBS covers 0.22 deg 2 of the COSMOS field 
and here we compare the z p hot determined between each 
survey. The NMBS catalog has similar data to the Ul- 
traVISTA catalog. It is also a K s -selected catalog and 
contains photometry from the Subaru broadband and 
medium band data. It also uses the same source-fitting 
code for the GALEX and IRAC photometry and the 
z p hot for the NMBS have been computed using EAZY. 
The primary differences are that the NMBS has 6 NIR fil- 
ters compared to 4 from UltraVISTA. This NMBS is also 
~ 0.6 mag shallower than UltraVISTA. The NMBS does 
have deeper optical photometry in the u*gVi'z' bands 
from the CFHTLS data which are not included in the 
UltraVISTA catalog. 

In the left panel of Figure 6 we plot the z p hot,uviSTA 
vs. z p hot, nmbs- Dark circles are for galaxies that are 
above the 95% mass-co mpleteness limit of the survey 
at a given redshift (see iMuzzin et al.l l2013f) , and light 
symbols are for those below the 95% mass-completeness 
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TABLE 4 

Summary of Photometric Redshift Catalog 



Column 


Parameter Name 


Description 


1 


id 


Object identifier number 


2 


z_spec 


Spectroscopic Redshift from zCOSMOS (no redshift = -1) 


3 


chi 


X 2 of the best-fitting template 


4 


z_peak 


Photometric redshift from the peak of the P(z) distribution 


5, 6 


168, u68 


Upper and lower 68% confidence intervals on z_peak 


7, 8 


195, u95 


Upper and lower 95% confidence intervals on z_peak 


9, 10 


199, u99 


Upper and lower 99% confidence intervals on z_peak 


11 


pcak_prob 


Peak probability 


12 


rrfit 


Number of filters used to determine z_peak 



2.0 



1.5 



Number of spectro: 5119 
Fraction >3a outliers: 80 (1.56%) 
.dz/(1 + z) = 0.013 




• Onodera et ol. (2012) 

• Bezanson et al. (2013) 
ovan de Sonde et al. (2013) 
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Fig. 5. — Left panel: Photometric redshift from the UltraVISTA K s -selected catalog vs. spectroscopic redshift from the zCOSMOS 
catalog (black dots). Only galaxies with high-quality spectroscopic redshifts and uncontaminated photometry are shown. Red galaxies with 
high-quality spectrosc opic redshifts from Onodera et a l. (2012) are plotted as red symbols. Also plotted are red galaxies with redshifts from 
Bezanson et al. (2013) and van dc Sande et al. (2012). Right panel: Difference between spectroscopic redshift and photometric redshift as 
a function of spectroscopic redshift. The dashed lines show the rms dispersion of 8z = 0.13/(1 + z), which is calculated from the zCOSMOS 
sample once catastrophic outliers are removed. The Onodera, Bezanson, and van de Sande spectroscopic redshifts suggest that the rms 
dispersion at z > 1.4 is likely to be larger than that measured from the zCOSMOS sample. 



limit. For the mass-complete sample, the agreement be- 
tween UltraVISTA and NMBS is excellent. The fraction 
of galaxies at < z < 4 that are > 5tr outliers is 2.0%, 
with a scatter in Sz/(l+z) = 0.026. Including the full 
sample of galaxies the number of > 5a outliers increases 
to 3.7%. 

In the right panel of Figure 6 we plot the difference 
between the z p hot as a function of the z p hoto,NMBS- At 
z > 1, there is a slight systematic bias between the two 
surveys of Sz/(l+z) = 0.005, with the z phot from NMBS 
being slightly higher for a given galaxy. This was also 
suggested by the small number of galaxies in this red- 
shift range with spectroscopic redshifts (§ 4.2). Small 
systcmatics like this are not completely unexpected. The 
NMBS used NIR medium bands in order to better-trace 
the Balmer/4000A bre ak and provide improv ed z p hoto for 
galaxies at z > 1.5 (sec Whitaker et al. 201 lj) . Even with 
the slightly deeper data available from UltraVISTA, it 
would be unrealistic to expect NIR broad bands to per- 
form as well as NIR medium bands when determining 

Z p hot • 

4.4. Correction for Extreme Blue Galaxies 

In the left panel of Figure 6, it is clear that a subsam- 
ple of the fainter galaxies have z p hot that do not agree 



well between the surveys. These are seen as the line of 
galaxies that have z p hot < 0.5 in UltraVISTA but a range 
of much higher z p hot in NMBS. There is a complemen- 
tary population running along the Y-axis; however, it is 
substantially smaller suggesting that on average the z p hot 
of this population may have been underestimated in the 
UltraVISTA catalog. 

As an example, we show the SEDs and P(z) of one of 
these objects in both surveys in Figure 7. Examination 
of the P(z) for these galaxies in both surveys shows that 
they are frequently bi-modal with a peak both at z < 
0.5, as well as a peak at z > 1.5. The SEDs of these 
galaxies are also typically very blue. They have a single 
feature, which is a break starting in the bluest optical 
bands. This break is interpreted as the Balmer-break in 
the low-redshift solution, and as the Lyman break for the 
high-redshift solution. 

These galaxies are problematic, but very few of them 
have spectroscopic redshifts so determining which solu- 
tion is the correct one is non-trivial. Several lines of 
evidence point to the NMBS solution of high-redshift Ly- 
man breaks being the correct ones. Firstly, many of the 
galaxies with the high-redshift solution do show weak 
Balmer breaks in the NIR medium bands (see Figure 7). 
These breaks are well traced by the NIR medium bands. 
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6. — Left panel: Photometric rcdshifts from UltraVISTA vs. photometric redshifts from the NMBS in the overlap region between 
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the two surveys. Galaxies with stellar masses > 95% mass-completeness limit for the survey are shown as dark circles, and those below 
the limit are shown as light circles. Right panel: Difference in photometric redshift between UltraVISTA and NMBS as a function of 
photometric redshift. The solid curve shows the median difference and the gray curves show the rms. The agreement between the two 
surveys is extremely good, particularly for the mass-complete sample where the >5 <r outlier fraction is 2.0% and the rms is <5z/(l+z) = 
0.026. 

This makes the u* band breaks clearer, and in many cases 
rules out the possibility of the break being the Balmer 
break, which is a weaker feature. Lastly, once rest-frame 
colors are computed (see § 5.2), it is clear that many of 
these galaxies lie in a portion of the UVJ diagram not 
populated by other galaxies. 

Given these lines of evidence, it is reasonable to expect 
that most of these galaxies are likely to be high-redshift 
galaxies and so we apply a correction to their photomet- 
ric redshift. The NMBS only covers a fraction of the 
UltraVISTA field, so we cannot use those data to im- 
prove things. Instead, we identify the population based 
on their location in the UVJ diagram. Galaxies that are 
much bluer than the overall z < 0.5 population in both U 
- V and V - J (rest-frame) are selected as the candidate 
high- z population. For these galaxies we re-run EAZY, 
but only allowing solutions in the range 1.5 < z < 6.0. 
This prevents a solution at z < 0.5 and forces the z p hot 
to the z > 1.5 solution. 

While this correction is somewhat ad hoc we note that 
it affects only a small fraction of galaxies in the catalog. 
Of the total sample of 262 615 sources, only 2415 (< 1%) 
are affected by this correction. More importantly, of this 
1%, the majority are quite faint, and as shown in Figure 
6, almost none lie above the mass-completeness limit of 
the survey no matter whether their solution is at high- 
rcdshift or low-redshift. Because of this, the correction 
of their z p hot does not affect any results based on the 
mass-complete sample. 

5. SED MODELING AND STELLAR POPULATION 
PARAMETERS 

5.1. Stellar Masses, Ages, and Dust Extinction 

Stellar population parameters are deter mined by fit- 
ting galaxy SEDs using the FAST code (|Kriek et al.l 
2009). We provide two catalogs of population parame- 
ters for the catalogs, one fit to SEDs generated using the 
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FlG. 7. — Example of a an extreme blue galaxy with a photo- 
metric redshift that is very different between the UltraVISTA and 
NMBS. The NIR medium bands in NMBS suggest a Balmer break, 
which is less well-resolved with the UltraVISTA broad bands. We 
have corrected the photometric redshift of the extreme blue galax- 
ies at z < 0.5 to their best fit solution at z > 1.5. The correction 
is primarily aesthetic, as these galaxies make up < 1% of all galax- 
ies and are typically well below the mass-completeness limit of the 
survey (at either redshift). 

They are also seen in the UltraVISTA broad bands, but 
are not as well traced, given that there arc less filters. 
As in Figure 7, the break in UltraVISTA is usually man- 
ifest as a single deviant point (either the Y or J band) 
to an otherwise good fit. Secondly, the NMBS uses the 
deep CFHTLS u* data which covers only a portion of the 
UltraVISTA field. It also uses as smaller color aperture 
(1.5" compared to 2.1"), so the photometric errors in the 
u* band in the NMBS are smaller by a factor of ~ 1.5. 
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TABLE 5 

Summary of Stellar Population Parameter Catalog 



Column 


Parameter Name 


Description 


1 


id 


Object identifier number 


2, 3, 4 


z, 168.Z, u68.z 


Photometric redshift and 68% confidence intervals from EAzy 


5, 6, 7 


ltau, 168Jtau, u68_ltau 


Best-fit value of log(r) and 68% confidence intervals 


8, 9, 10 


lage, 168Jage, u68_lage 


Best-fit value of log(t) and 68% confidence intervals 


11, 12, 13 


A„, 168_A„, u68_A„ 


Best-fit value of A„ and 68% confidence intervals 


14, 15, 16 


lmass, 168_lmass, u68_lmass 


Best-fit value of log(M s t ar /M0) and 68% confidence intervals 


17, 18, 19 


lsfr, 168Jsfr, u68_lsfr 


Best-fit value of log(sfr) from the SED and 68% confidence intervals 


20, 21, 22 


lssfr, 168_lssfr, u68_lssfr 


Best-fit value of log(ssfr) from the SED and 68% confidence intervals 


23 


chi 


X 2 of best-fitting model 
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Fig. 8. — Comparison of stellar masses determined in the NMBS 
and UltraVISTA as a function of stellar mass. The dark solid curve 
is the running median and the lighter solid curves encompass the 
68-percentile of the distribution There is a trend for systematically 
higher masses in the NMBS, but no significant trend with stellar 
mass. 



iBruzual fc Charlotl (2003) mod els an d one fit to SEDs 
generated using the iMarastonl (]2005l ) models. For both 
sets o f models we assume solar metalli city, a IChabrierl 
(|2003h IMF, and a lCalzetti efaLl ()2000f ) dust extinction 
law. 

To construct the set of template SEDs, we assume 
galaxies have exponentially-declining star formation his- 
tories (SFHs) of the form SFR oc exp(-i/r), where t is 
the time since the onset of star formation and r is the 
e-folding star formation timescale in units of Gyr. We 
allow Log(r) to vary between 7.0 and 10.0 in increments 
of 0.2, and Log(i) to vary between 7.0 and 10.1, in in- 
crements of 0.1. For all galaxies we restrict t to be less 
than the age of the universe at the redshift of the galaxy. 
We also fit for visual attenuation of the galaxies (A^) as- 
suming a uniform dust screen geometry and allow A„ to 
vary between and 4. All galaxies are fit assuming their 
redshift is the best-fit EAZY z p h ot . In all, we fit four pa- 
rameters per galaxy: r, t, A„, and a normalization. The 
stellar mass (M star ) is then determined from mass-to- 
light ratio of the best-fit SED multiplied by the best-fit 
normalization of the SED. The layout of the catalogs of 
best-fit stellar population parameters is summarised in 
Table 5. 

In Figure 8 we show a comparison between the M star 



measured for the galaxies in the NMBS catalog vs. the 
M star measured for galaxies in the UltraVISTA catalog 
as a function of mass. In general the agreement is good. 
There is a systematic difference of 0.05 dex between the 
two with galaxies being more massive in the NMBS cat- 
alog. 

5.2. Rest-frame Colors 

The rest-frame U - V vs. V - J diagram has 
become a popular way to differen tiate between star- 
forming and quiescent galaxies (lLabbe et al.l 120051: 
Wuvts et al.l 120071: IWiTliams et al.ll2009l: iBrammer et al l 



20111: IWhitaker et all 12011b iPatel et"afl l2012al lbl). This 



approach is similar to the obs erved BzK diagram that 
has been used in the past (e.g-. lDaddi et al.l 120051 120071 
and numerous others) , but allows for a cleaner separation 
of star-forming and quiescent galaxies because colors are 
defined in the rest-frame rather than observed frame, re- 
moving any redshift-dependence of the colors. 

We calculate rest-frame U - V and V - J colors for all 
galaxies using EAZY. EAZY determines the colors by 
integrating the best-fit SED through the redshiftcd fil- 
ter curves over the appropriate wavelength range. For 
the U and V fil t er we use the response curves defined in 
iMai'z Apellanizl (|2006|) . an d for the J filte r we u sed the 
2MASS filter curve from ISkrutskie et al.l (|2006l ). The 
rest-frame U - V and V - J colors are listed in Table 6. 



5.3. The UVJ Diagram 

In Figure 9 we plot grayscale histograms of the U - V 
vs. V - J colors for galaxies with S/N(K S ) > 7 in various 
redshift bins between < z < 3.5. The galaxy popula- 
tion is clearly separated into two clumps in color-color 
space up to z = 2. The reddest of the two clumps has 
colors similar to that of quiescent and passively-evolving 
galaxies, whereas the bluer clump has a much wider range 
of colors, usually interpreted as the star-forming popula- 
tion with a range of dust extinctions and geometries (e.g., 
Labbe et al.l 120051; IWilliams et al.l |21M IBrammer et al.l 



2009HPatel et al.ll2012bh . In general, the UVJ diagram in 



COSMOS shows t he same structure as has been seen in 
other surveys (e.g.. IWilliams et aI1l2009t [Whitakc r et al.l 
|2011| ): however, the superior volume in UltraVISTA pro- 
vides an increase in the number of galaxies by a factor 
of- 4- 5. 

An interesting feature of the diagram revealed by the 
large sample of galaxies in UltraVISTA is the continued 
reddening of the quiescent population with decreasing 
redshift. Whereas the quiescent population is primarily 
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Fig. 9. — The U - V vs. V - J color-color diagram for galaxies with S/N(K S ) > 7. The well-known bi-modality between quiescent galaxies 
and star forming-galaxies can be clearly seen in the galaxy population up to z ~ 2. Thereafter the bi-modality becomes less pronounced. 
Also visible in the figure is the continued reddening of the passive population with decreasing rcdshift. This reddening is expected from a 
passively-evolving population. 



TABLE 6 

Summary of Rest-frame color and UV and IR Luminosity Catalog 



Column 


Parameter Name 


Description 


1 


id 


Object identifier number 


2 


U - V 


Rest-frame U - V color 


3 


V - J 


Rest-frame V - J color 


4, 5, 6 


L2800 


Total UV luminosity and 68% confidence intervals 


7, 8, 9 


L/ii 


Total IR luminosity and 68% confidence intervals 


10, 11, 12 


SFRfJV, uncorr 


UV star formation rate and 68% confidence intervals, uncorrected for dust extinction 


13, 14, 14 


SFR IR 


IR star formation rate and 68% confidence intervals 



located in a single clump with colors of U - V — 1.7 and 
V - J = 1.0 at z = 2, it appears as more of a sequence 
at z = 0. This reddening with decreasing redshift is 
precisely what is expect ed for a passively- evolving pop- 
ulation of galaxies (e.g.. lWuvts et al.ll2007| ). 

It is tempting to make further inferences on the evolu- 
tion of the galaxy population using Figure 9. Doing so 
requires a more quantitative analysis, in particular it is 
important to adopt the appropriate limits in M star with 
rcdshift in order to define a mass-complete sample. Fig- 
ure 9 is meant to be illustrative and is made with a S/N 
cut, not a M stQr -cut, so it becomes increasingly skewed 
to lower-mass galaxies at lower- redshift. A full analysis 
of the UVJ diagram using the appropriate mass limits 
will be presented in a future paper. 

5.4. Example SEDs 



As an example of the quality of the SEDs and SED fit- 
ting, we show examples of the observed galaxy SEDs and 
the best-fit EAZY SEDs from the catalog in Figures 10 
-12. At each redshift, SEDs classified by their location 
in the UVJ diagram are shown in order to demonstrate 
the SEDs for a range of different galaxy types. We chose 
a total of five different SED classifications, two types of 
SEDs in the quiescent population, and three types in the 
star-forming population. Within the quiescent popula- 
tion we chose those that are located near the red-tip in 
U - V vs. V - J (labelled "old"), and those that are near 
the blue-tip in U - V vs. V - J, but still have quiescent 
colors. These blue-tip galaxies have colors that are simi- 
lar to recently quenched galaxies, or "post star formation 
galaxies" fe.g- lKriek et^[2Tjia IWhitaker et al.ll2012aft 
and are labelled "post-SF" . 

For the star forming population, we divided the se- 
quence of galaxies going from blue U - V and blue V - 
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Fig. 10. — Example SEDs for galaxies of different types at < z < 1.0 with S/N(K S ) > 10 chosen from their location in the UVJ diagram 
(see text). In general the fits are very good, and clear features such as the Lyman break, Balmer break, and 4000A break can be seen in 
the galaxy populations. 



J colors, to red U - V and red V - J colors into three 
bins. The three regions cover the bluest population, an 
intermediate population, and the reddest (dustiest) pop- 
ulation (labeled blue-SF, inter-SF, and dusty-SF, respec- 
tively). 

The galaxies in each bin are chosen to have S/N(K S ) 
> 10, but are otherwise selected at random so that they 
are representative of typical SEDs for that SED type and 
redshift range. The only exception to this is that for the 
lowest redshift bin we use only galaxies that have z spec 



from zCOSMOS in order to show the agreement between 
photometric and spectroscopic SEDs. 

Figures 10-12 demonstrate the excellent quality of the 
SEDs and SED fitting, all the way up to z = 2.5. They 
also show how at z < 1 the optical medium bands trace 
the Balmer and 4000A break extremely well, and why 
such good quality z p hot can be determined. At z > 1.5 
the UltraVISTA bands trace the break. In particular, 
it can be seen how the Y-band data is extremely use- 
ful for connecting what would be a large gap in wave- 
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Fig. 11. — As Figure 10, but for galaxies at 1 < z < 2. 



length space between the z'-band and J-band. At z > 
2.5 the IRAC data becomes increasingly important, be- 
cause these are the only filters that remain redward of 
the break. The best-fit SEDs for all galaxies using EAZY 
and and FAST for both BC03 and M05 models are avail- 
able for download with the catalog. 

5.5. UV and Infrared Luminosities and SFRs 

The SED fitting to the BC03 and M05 models re- 
sults in an estimated SFR for each galaxy (see Table 
5). Those values can be considered total, dust-corrected 



SFRs (primarily constrained by the UV flux) . They may 
be indicative; however, we caution against using those 
values in a quantitative way. This is because the esti- 
mated SFR is strongly influenced by the assumption of 
an exponentially-declining SFH. If this is not the cor- 
rect SFH (which i s likely to be true for many star form- 
ing g alaxies, e.g., iMaraston et al.l l2010t iPapovich et al.l 
120111 ). then they can be substantially in error. 

A better measure of the instantaneous SFR comes from 
the rest-frame UV flux (L2800) and the rest-frame total 
infrared luminosity (Lir). With the availability of rest- 
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Fig. 12. — As Figure 10, but for galaxies at 2 < z < 3. 

frame UV information (from GALEX to optical) as well 
as rest-frame MIR data (from MIPS-24/im) we can cal- 
culate these quantities, or upper limits on them, for all 
galaxies. These can then be converte d to SFRs based on 
the standard conversion factors fe.g.. lKennicuttj [l998'). 

To determine L2800 we use EAZY. Similar to the rest- 
frame colors (§5.2), EAZY integrates the best-fit tem- 
plate over the wavelength rang e 2600A - 2950A to deter- 
mine an L2800 for all galaxies. Because there is GALEX 
data, this quantity is constrained by data over the full 
rcdshift range of the galaxy population. The L2800 values 



are listed in Table 6. 

To determine the Lir, we ex trapolate the m e asure d 
24/xm flux usi n g tem plates. The lCharv fc Elbad (|2001[ k 
iDale fc Heloul |2002t ) models are based on local galaxy 
templates calibrated using IRAS and implicitly assume a 
correlation between Ljr and the dust temperature (T^). 
As has been discussed in recent papers, it appears that 
this c orrelation may not hold up to the high est redshifts 
(e.g.. iMuzzin et al.l 120101: lElbaz et al.ll2010t ). Instead of 
using the luminosity-dependent templates, we use a sin- 
gle template to determine Ljr for all galaxies. This ap- 
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Fig. 13. — Log(SFR) vs. M s tar in difference redshift bins for star forming galaxies selected using the UVJ diagram (Figure 9). As in 
many previous works th ere is a corre l ation betw een the SFR of gala xies and their mass. The blue and red curves are power-law fits to the 
relations determined by Elbaz ct al. (2007) and Daddi ct al. (2007) using data from the GOODS survey. The normalization and slope of 
the star formation main sequence evolves similarly between the two studies. 



proach has been advocated by many recent studies (e.) 
IWuvts et al.ll2008l : lMuzzin et alfe OlO; El baz et alJ^Olt 
We use the log-average of the [Dale & Hclou ( 2002) tem- 
plates for this computation (see, IWuvts et al.ll2008[). but 
note t hat using the log-average of the ICharv fc Elbazl 
20011) templates w ould provide very similar results 



Muzzin et al.l [2010h . We list the Lm determined us- 
ing this method in Table 6. 

We convert the L2800 into a SFKuv.uncorr using the 
conversion factor SFRjjy uncorr = 3.234 x 10 -10 L2 



^2800 



from IKennicutt (|1998). adapted to a Kroupa IMF by 
iBell et al.l(|2005l ). We note that this is the observed SFR, 
and is not corrected for dust-extinction. The Fj R is con- 
verte d into a SFR rR u sing SFR IR = 0.98 x 10~ 10 L IR 
from IKennicutt (|1998l) . adapted to a Kroupa IMF by 
IBell et al.l (|2005l ). The total SFR of the galaxy can then 
be determined via SFRtot = SFKuv,uncorr + SFKir. 
These values are also listed in Table 6. 



5.6. The Star Formation Main Sequence 

Many recent studies have shown evidence for a cor- 
relation between the SFR of star forming galaxies and 
their M star (e.g.. iNoeske et al.l 120071: lElbaz et al.l 120071: 
Daddi et al.l l2007i IWuvts et alj 1201 It IWhitaker et al l 
2012bD . This correlation has become known as the "star 
formation main sequence" . In Figure 13 we plot greyscale 
histograms of the star formation main sequence for star 
forming galaxies in the UltraVISTA catalog between 0.0 
< z < 2.5. Star forming galaxies have been selected by 
their loca tion in the UVJ diagr am using the prescriptions 
defined in Mu zzin et al.l (|2013l ) , which are based on those 



determined bv lWilliams et al.l (|2009l ). The SFRs plotted 
in Figure 13 are total SFRs which have been computed 
from SFR tot = SFR[/v>„ corr + SFR//;. Only sources 
that have a > 3cr detection at 24/Ltm are shown. 

Figure 13 shows clearly the existence of a main se- 
quence of star formation for galaxies in the UltraVISTA 
catalog. As a comparison we plot the power-law fits to 
the main se quenc e at z ~ 1 and z ~ 2 measured by 
lElbaz et al.l (|2007l) and lDaddi et al.l (|2007l) . respectively 
using the GOODS data. There is reasonable agreement 
between the evolution of the normalization of the rela- 
tions between the two datasets. There may be slightly 
different slopes and normalizations, but this may be be- 
cause the redshift ra nges shown in Figu re 13 are not exact 
match es to those in lElbaz et al.l (|2007l ) and lDaddi et al.1 
((20071) . A more detailed look at the star formation main 
sequence from UltraVISTA will be presented in a future 
paper. 

6. CONCLUSION 

In this paper we have presented a public K s -selected 
catalog covering 1.62 deg 2 of the COSMOS/UltraVISTA 
field. The photometric catalog consists of PSF-matched 
photometry in 30 bands and covers the wavelength range 
0.15/^m — s> 24/^m. The catalog contains 262 615 sources 
down to the 3cr limit of K s (2.1") < 24.35, and 179 291 
sources down to the 90% completeness limit K s to t < 
23.4. 

Photometric redshifts have been computed for all 
galaxies using the EAZY photometric redshift code. 
Comparison of the z p hot with z spec from zCOSMOS and 
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other spectroscopic surveys shows that at z < 1.5 the 
catastrophic outlier fraction is low (1.56%), and the rms 
scatter is low (Sz/(l+z) = 0.013. The z p hot also agree 
well with Zphot determined using the NMBS. 

Included with the catalog is a set of M star and stellar 
population parameters computed using the FAST SED 
fitting code. These M sta r show good agreement with 
those calculated in the NMBS, with only a 0.05 dex sys- 
tematic difference. Rest-frame U - V and V - J colors 
are computed for all galaxies. The population of galax- 
ies in the COSMOS/UltraVISTA field shows a clear bi- 
modality in color-color space up to z ~ 2, and thereafter 
the bi-modality begins to break down. 

The catalog contains measures of L2800 and L/^j as well 
the inferred SFRs from those parameters. Plotting these 
against the M star for UVJ-selected star forming galaxies 
shows that there is a star formation main sequence up to 
z ~ 2.5. The evolution of the main sequence is consistent 
with previous measurements. 

In a companion paper to this one uzzm "eT^l2013h 
we show the evolution of the stellar mass functions of 

11 http:/ /www. strw.leidenuniv.nl/galaxyevolution/ULTRAVISTA/ 



galaxies up to z ~ 4 using the K s -selected UltraVISTA 
catalog. The photometry, photometric redshifts, stel- 
lar population parameters, rest-frame colors, UV and IR 
SFRs from the K s -selected catalog are now made avail- 
able to the astronomical community though the catalog 
website 11 . We hope that it proves to be a useful resource 
for further galaxy evolution studies. 
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APPENDIX 

COMPARISON OF ULTRAVISTA K s ZEROPOINT TO OTHER DATASETS 

In this appendix we make a comparison between the K s -band photometry between several datasets that cover the 
COSMOS field to test the photometric zeropoint of the UltraVISTA K s imaging. For the catalog, the zeropoint in 
K s is important because K s is used as the anchor filter when calculating the zeropoint offsets in the remainin g filters 
(§ 4.1). We compare the UltraVISTA KVband to the K-band imaging from the NMBS (IWhitaker et al]|2011h . which 
reaches a 5cr depth of K = 23.5. We also compare the CFHT / WIRC AM K s imaging from iMcCracken et al] (|2010ft . 
which in th e deepest regions reach es a 5a depth of K s = 23.65 in a 2" aperture. Lastly, we compare to sources detected 
in 2MASS (|Skrutskie et al]|2006l ). although we note that the overlap in in dynamic range between UltraVISTA and 
2MASS is limited. 

For comparison with the WIRC AM and NMBS data we use the region of the survey that is in common to all 
three datasets. This region covers ~ 0.22 deg 2 and is located in the northwest corner of the UltraVISTA field. We 
perform PSF matching on all three datasets, matching to the worst-seeing image, which is the NMBS (~ 1.1"). Ten 
bright unsaturated stars in the region are used in order to determine the convolution kernel. Once the PSFs are 
homogenized we use SExtractor to find objects on each image individually and measure fluxes within a 2.1" diameter 
aperture. Object detection is performed on each image separately because there are small, but noticable differences 
in the astrometry between the surveys. If object detection was performed with SExtractor in dual-image mode, small 
astrometry differences could appear as zeropoint differences. Source lists from each catalog are matched using a 0.5" 
search radius. Using a small radius means that some matches may be missed, but the tradeoff is that there are few 
ambiguous matches. 

In Figure 14 we plot the comparison of the UltraVISTA K s photometry vs. the WIRC AM and NMBS K s photometry. 
In general, the comparison is good. There is a small systematic difference between the surveys; however, there is no 
systematic difference as a function of magnitude. If we compare the differences in the magnitude range 16 < K s < 20, 
we find an offset between UltraVISTA and WIRCAM of 0.051 mag, and an offset between UltraVISTA and NMBS of 
0.078 mag, with UltraVISTA being fainter in both cases. 

The difference o f 0.05 4 mag between the UltraVISTA and the WIRCAM data is very similar to that measured by 
McCrack en et al.l (|2012f) (see their Figure 10) for BzK-selected stars in the fields of the two surveys. The difference 
between UltraVISTA and the NMBS is slightly larger than the difference between UltraVISTA and WIRCAM. That 
difference implies that the WIRCAM photometry should be 0.027 mag f ainter than the NMBS photometry. Indeed, 
this implied difference is very similar to the actual difference measured by IWhitaker et aTl (|201 lh who performed PSF 
matched photometry between the NMBS and the WIRCAM data (see their Figure 12). 

In the bottom panel of Figure 14 we also compare the UltraVISTA photometry to photometry from the 2MASS 
catalog. For this comparison we do not perform PSF matching, but match objects between the surveys over the full 
1.62 deg 2 UltraVISTA field. We compare the total magnitudes measured for UltraVISTA (before any zeropoint offset 
has been applied), to the "default" magnitudes of objects extracted from the 2MASS point source catalog. As Figure 14 
shows, there is reasonable agreement; however, the dynamic range over which the comparison can be made is limited. 
If we compare objects at 15.5 < K s < 17.0 we find a systematic difference of 0.054 mag, with again UltraVISTA being 
slightly fainter. 

Although it is not completely clear which zeropoint is the most trustworthy, these comparisons do suggest that the 
DR1 UltraVISTA K s band zeropoint may be too large at the 0.05 - 0.08 mag level. In order to provide consistency 
with our previous work, we have elected to adjust the zeropoint, making it 0.08 mag brighter so that it matches the 
NMBS zeropoint. We note that this offset in zeropoint does not have a significant impact on the colors in the catalog 
because the zeropoint offsets in the other bands are derived relative to the K s -band. This means that parameters 
derived from the colors such as the z p hot are also unchanged. The best-fit M/L ratios are also unchanged but the 
scaling of the K s band means that the M star are systematically changed by adjusting the K s zeropoint; however, this 
change is quite small. The offset of 0.08 mag corresponds to a difference of 0.03 dex in M stQr . If we were to apply 
the zeropoint implied by 2MASS or the WIRCAM data, it would imply stellar masses that are only 0.01 dex different 
than those in the current catalog. 
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Fig. 14. — Top two panels: Comparison between K s photometry in UltraVISTA to the K s and K photometry from WIRCAM and the 
NMBS. Magnitudes have been measured in a 2.1" aperture on PSF matched images. Bottom Panel: Comparison of the UltraVISTA K s 
photometry to 2MASS K s photometry over the full UltraVISTA area. Blue dotted lines indicate the magnitude range used to compute the 
median offset, and the red solid line shows the median offset. There is a < 0.1 mag difference between UltraVISTA and the other surveys; 
however, UltraVISTA is systematically too faint and so we have elected to redefine the K s zeropoint to match the NMBS zeropoint. 



